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Abstract 

Ultrafast laser-induced demagnetization of Gd(OOOl) has been investigated by magneto-induced 
optical second harmonic generation and the magneto-optical Kerr effect which facilitate a compar- 
ison of surface and bulk dynamics. We observe pronounced differences in the transient changes of 
the surface and bulk sensitive magneto-optical signals which we attribute to transfer of optically 
excited, spin-polarized carriers between surface and bulk states of the Gd(OOOl) film. A fiuence 
dependent analysis of the bulk magnetization dynamics results in a weak variation of the demagne- 
tization time constant, which starts at about 700 fs and increases by 10% within a fiuence variation 
up to 1 mJ/cm^. We compare these results with fiuence dependent changes in the transient energy 
density calculated by the two temperature model. The determined characteristic times of excess 
energy transfer from the electron system to the lattice, which is mediated by e-ph scattering, range 
from 0.2 - 0.6 ps. Such a more pronounced fiuence dependent change in the characteristic time 
compared to the observed rather weakly varying demagnetization times suggests a more advanced 
description of the optically excited state than by the two-temperature model. 

PACS numbers: 78.47.J-, 71.38.-k, 75.70.Ak, 78.70.Dni 
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I. INTRODUCTION 



The understanding of spin dynamics in ferromagnets is widely based on radio frequency 
spectroscopy^ and inelastic scattering techniques employing photons^ or neutrons^. These 
experimental methods provide frequencies and decay rates of magnetic excitations by the 
line position and line width, respectively. More recently time domain techniques that em- 
ploy short laser pulses have been established in a considerable number of laboratories and 
provide a complementary approach to analyze spin dynamics^"-. In pump-probe experi- 
ments the magnetization dynamics is driven by an external stimulus like a laser or magnetic 
field pulse and is probed by a second, time-delayed laser pulse. Using femtosecond laser 
pulses, which are nowadays routinely provided by Ti: Sapphire laser systems, the field of 
femtosecond magnetization dynamics has developed enormously. One reason for the inter- 
est in this approach might be that it is not only complementary to the established tools for 
magnetization dynamics, but that the non-equilibrium state and the respective dynamical 
magnetization processes have become accessible for systematic experiments^. 

A number of novel phenomena like laser-induced demagnetization on sub-picosecond time 
scalea^i^, laser-driven formation of ferromagnetic order-, and magnetization reversal by indi- 
vidual laser pulses^ demonstrates the potential for spin manipulation in solid state materials 
on ultrafast time scales. At the same time, these phenomena challenge our understanding 
of ferromagnetism and ask for an appropriate description of the observed dynamics. 

Beside the investigation of multi-constituent materials with a potential for novel phenom- 
ena, it is essential to develop insight into the more general laser-driven dynamics to establish 
this understanding. Finally, the meanwhile well established sub-picosecond demagnetization 
of elemental metallic ferromagnets is still controversial because the underlying elementary 
interactions remain under discussion. It has been shown recently that descriptions based 
on phonon-mediated spin-flip scattering (Elliott- Yafet tjpe)-^^, electron-spin-flip scatter- 
ing (Stoner excitation)^^, and ballistic transport of spin-polarized charge carriers^^ can be 
used to describe rather similar experimental data of femtosecond demagnetization in Ni 
films. In the phonon-mediated process the magnetic moment is considered to be transferred 
directly to the lattice through phonons with appropriate symmetry obeying angular momen- 
tum conservation. In the electronic spin-flip scattering process the magnetic moment is not 
transferred within the elementary scattering process itself, but is dissipated by secondary 
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low energy spin excitations, which in their sum lead to demagnetization^^. In the transport 
concept^ the spin polarization is basically taken out of the ferromagnet into a substrate or, 
alternatively, out of the experimental observation window. Therefore, the current under- 
standing is not conclusive and a discrimination of these above processes is important at this 
stage. 

Our efforts have concentrated on rare earth ferromagnets and in particular on Gadolinium 
which is considered to be a model system for a Heisenberg ferromagnet with a magnetic 
moment localized at the ion core. We have employed in earlier studies time-resolved surface 
sensitive techniques which are optical magneto-induced second harmonic generation (SHG) 
and photoelectron spectroscopy, to analyze the excited state and its relaxation. The response 
to the laser excitation is characterized by a coherent phonon-magnon mode localized at the 
surface and incoherent dynamics of optically excited electron-hole pairs, lattice vibrations, 
and magnetic excitations. A comprehensive overview is given in Refs.— >i^. 

Here we present experimental results obtained by the femtosecond time-resolved magneto- 
optical Kerr effect (MOKE) on Gd(OOOl) films to analyze the magnetization dynamics in the 
bulk part of the Gd film, which we compare with the surface sensitive SHG signal. We find 
pronounced differences during the first 2 ps which we attribute to transfer of spin-polarized 
carriers between the surface and bulk part of the film. Moreover, we study the fluence 
dependence of the demagnetization process and compare it with predictions based on the 
two-temperature model. 

II. EXPERIMENTAL ASPECTS 

The experiments were carried out at a setup that combines an ultrahigh vacuum (UHV) 
chamber and a femtosecond (fs) Ti: Sapphire laser. For details se©^. Gd was evaporated 
from an electron beam heated crucible under UHV conditions onto a W(llO) substrate at 
300 K. Annealing to 700 K subsequent to the deposition results in smooth epitaxial filmsi^. 
Here we studied films with 20 nm thickness. Pump-probe experiments were performed at 
50 K equilibrium temperature of the sample which was cooled by a liquid He cryostat. Laser 
pulses of 35 fs duration and 40 nJ energy were generated by a cavity dumped Ti:Sapphire 
oscillator at 800 nm central wave length with a repetition rate of 1.52 MHz. For pump-probe 
experiments the laser output was divided at a ratio 4:1. To detect the p-polarized second 
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harmonic (SH) generated by the p-polarized probe pulse a monochromator selecting 400 nm 
and a photomultiplier were used. To analyze the MOKE a balance detection scheme was 
used to measure the polarization rotation of the probe pulse 9 in the longitudinal MOKE 
geometry. The magneto-optical signals were analyzed by a lock-in amplifier for MOKE and 
a single photon counter for SHG to detect the differential signals for open and blocked pump 
in opposite saturation magnetic fields {Hl'^) as a function of pump-probe delay t. Formally 
negative delays t < denote the state before the optical excitation has occurred. For MOKE 
we show below 

Oo Oo ^{hI, t < 0) - ^{H^, t < 0) 

where i) is the angle between the linear polarization of the refiected probe pulse and a fixed 
reference polarization, 6 is the magneto-optical Kerr rotation. is the time-dependent 
change in the Kerr rotation which represents the magnetization of the bulk part of the Gd 
film. 

SHG is detected in the transversal geometry leading to magneto-induced changes in the 
SH intensity. 

It « (EZn" + (EfJ) ' ± COS 0. (2) 

Here E^^^ and -£^2^^^ SHG optical fields which behave as even or odd with respect 

to reversal of the magnetization M. The phase between these two field contributions is </> 
which is smaller than 15° and weakly t ime- dependent . Since E^^"^ oc M we plot below 

. odd.,^ ^ - < 0) ^ Mjt) - Mjt < 0) _ AM(t) 



{t < 0) M{t < 0) M, 







which represents the time-dependence of the surface sensitive magneto-optical signal. 



III. EXPERIMENTAL RESULTS AND DISCUSSION 

The article presents results of time-resolved SHG and MOKE. The first section focuses 
on the comparison of surface and bulk sensitive detection of the transient magnetic state in 
Gd(OOOl). The second section contains a detailed study of time-resolved MOKE and reports 
on the fiuence dependence of the bulk demagnetization. 
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A. Surface and bulk demagnetization dynamics 

We start by comparing the transient SHG and MOKE signals that probe the surface 
and bulk of the Gd(OOOl) film, respectively. In Fig. [H top panel, we show A^^^ (solid 
line, left axis, surface) and A6/6q (dotted line, right axis, bulk) taken at nominally identical 
experimental conditions. Both signals indicate a pronounced laser-induced demagnetization, 
but their transient behavior is very different. The magnetic SHG signal is reduced within the 
laser pulse duration of 35 fs by 0.5 and exhibits on that level oscillations that are damped 
out within 3 ps. The signal starts to recover after 1.5 ps and returns to its value before 
optical excitation after several 100 ps (not shown, see Ref.— ). The oscillatory part has 
been explained before as a coupled phonon-magnon mode localized to the surface^>^. It is 
therefore expected, that it is not found in the bulk signal. The MOKE data can be described 
by a continuous reduction that will be fitted by a single exponential decay below and tends 
to saturate at a delay of 3 ps, before it demagnetizes further, not shown see Ref. — . Overall, 
the bulk sensitive signal is reduced on a much slower time scale than the surface one and is 
in agreement with femtosecond x-ray magnetic dichroism studies^^ if the different excitation 
densities are considered. We therefore conclude that the time-resolved MOKE signal probes 
the time-dependent magnetization of the film. We note here that a detailed comparison of 
MOKE ellipticity and rotation points to time-dependent variations of the magneto-optical 
constants at delays where the excess energy resides predominantly in the electronic system, 
as will be discussed in a forthcoming publication^i. 

To explain the difference between the surface and bulk demagnetization it is informative 
to inspect the electronic structure at the surface and to take transfer processes into account. 
Fig. [H bottom, depicts several data sets which in their combination represent the electronic 
valence states of epitaxial Gd(OOOl) films. For a ferromagnetically ordered situation the 
states are exchange-split due to intra-atomic exchange interaction with the strong magnetic 
moment of the half filled 4/ shell. At the F-point the bulk 5(i-states appear at 1.4 eV (mi- 
nority) and 2.4 eV (majority) binding energy and disperse towards Ep with increasing k\^. 
The majority component of the exchange-split 5(i^2-surface state is dominantly occupied 
and the minority one unoccupied. The system also contains unoccupied exchange-split bulk 
states^^. 

For linear polarized light, within the dipole approximation, optical transitions proceed 
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FIG. 1: Top panel: Time-dependent magneto-optical signals measured on 20 nm thick Gd(OOOl) 
films, which were grown epitaxially on W(llO). The solid line shows the pump-induced change 
of the magnetic SHG contribution ^2^'^ sensitive to the surface. Similar results were published 
before^»2ii2^. The dotted line depicts the pump-induced MOKE polarization rotation AO normal- 
ized to the static rotation 6q. The bottom panel is reprinted from Ref.^^ and shows the valence 
electron states of Gd(OOOl) taken from normal direction photoemission (circles), inverse photoe- 
mission (solid lines), and scanning tunneling spectroscopy data above the Fermi level (solid line, 
filled)2^. The exchange-split surface state (filled area) appears around the Fermi level. Vertical 
arrows represent the majority and minority character of the electronic states. Indicated are the 
two main absorption channels for 1.5 eV pump photons and the resonant second harmonic probing 
scheme. 



among occupied and unoccupied electronic states in (bulk) valence bands with the same 
spin character. In case of the excitation by the 800 nm (or 1.55 eV) pump laser pulses, an 
additional excitation channel becomes available at the surface because of resonant transitions 
coupling the bulk and the localized surface electronic states, see Fig. [T], bottom. Therefore 
a difference of the surface sensitive SHG response with respect to the bulk sensitive MOKE 
can be expected. Similarly to a charge-transfer excitation across molecular interfacea^I, such 
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surface-bulk resonances lead to effective electron transfer between surface and bulk which 
modify the transient population of the minority and majority components of the surface 
state. The population of the initially occupied majority surface state component is reduced 
due to the surface-to-bulk spin- up electron transfer, see Fig. [H bottom. At the same time, the 
population of the initially unoccupied minority surface state component increases due to the 
bulk-to-surface spin-down electron transfer, which can alternatively be viewed as a surface- 
to-bulk spin-down hole transfer. The carriers in the surface state which are transfered to bulk 
bands represent initially a wave packet at the surface which spreads and propagates into the 
bulk material with its Fermi velocity vp of about 1 nm/fs^. Thus, a spin-polarized current 
which propagates from the surface to the bulk is optically excited. In the vicinity of the 
surface this current propagates ballistically and consists of electron and hole contributions. 
Spin-down holes have the same spin polarization as spin-up electrons. Therefore, the spin 
components of these two contributions will add up. If the charge of electrons and holes is in 
sum zero we conclude that a net spin current between surface and bulk is excited. It is likely 
that the efficiency of electronic transitions in spin-up and spin-down channels are different 
(Fig. [H bottom) which will lead to a prevalence of either the electron or hole component of 
the spin-polarized current. Note that this does not necessarily imply charging of the surface 
because a transient charge imbalance will be screened by electron rearrangement on the time 
scale of the inverse plasma frequency^. Note that on the basis of these considerations an 
effect due to optically excited transfer of spin polarization between surface and bulk should 
proceed within few femtoseconds, i.e. well within the time resolution of the SHG experiment. 

As seen in Fig. [H top panel, the pronounced initial drop in the surface sensitive magneto- 
optical signal occurs within the experimental time resolution and clearly faster than the bulk 
demagnetization. Therefore, we explain the pronounced difference in the bulk and surface 
demagnetization times by spin transfer between surface and bulk of the film. The time scale 
at which the magnetic SHG signal A^^^ changes from to -0.5 at time zero is in agreement 
with a ballistic character of these transport effects. 

The pump-induced reduction of the relative magneto-optical signals differ by a factor of 
six. We consider the resonant surface excitation among bulk and the surface states near F, 
see Fig. [H as essential for an explanation of the pronounced effect at the surface. 
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FIG. 2: Symbols indicate experimental data of time-resolved MOKE for different relative pump 
fluences F/ Fq , with Fq ~ 1 m J /cm^ . The solid lines are fits considering a single exponential decay. 

B. Fluence dependent demagnetization dynamics in the bulk 

The setup used for the present study allowed in addition a more detailed analysis of 
the bulk demagnetization phenomenon. Using a combination of half wave plate and Glan- 
Thomson polarizer the absorbed fluence has been reduced step wise from the maximum 
value F / Fq = 1, Fq was determined to be 1.0 ± 0.3 mJ/cm^. Fig. |2] shows representative 
time-resolved MOKE curves for different F/Fq up to delay times of 4 ps. The values at 
4 ps decrease linearly with F/Fq. This ensures that we are analyzing a low excitation den- 
sity regime reasonably far away from a full demagnetization of the sample, where magnetic 
fluctuations would contribute to the ultrafast magnetization dynamics^^'^. Near time zero 
the data feature a less clear behavior and an effectively positive contribution in AO/Oq. The 
time delay at which the signal crosses the zero line shifts with increasing fluence closer to 
time zero. While for F/Fq = 0.2 the amplitude of the positive and negative contributions 
are comparable to each other, the negative one dominates for larger F/Fq. In this study we 
focus on the pronounced demagnetization dynamics, i.e. at the fluences where the negative 
contribution dominates. The origin of the positive contributions are currently under discus- 
sion since they could originate from transient changes in the magneto-optical constants or 
the spin polarization of the 5d electrons. 

All data were fitted by a single exponential time dependence with fixed time zero and 
variable AO/Oq at and 4 ps. The obtained fits are shown in Fig. [2] and describe the 
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FIG. 3: Characteristic times determined from the decay times in the single exponential fitting of 
A9/9q as a function of relative fluence. The line is a linear fit to the data with F/Fq > 0.30 

experimental data well. Fig. |3] depicts the time constants determined by the fitting 
procedure. The numbers range from about 0.5 to 0.8 ps with a trend towards larger times 
for higher F/Fq. The smallest are obtained for the lowest F/Fq = 0.20 and 0.27 and 
the respective Tm are out of the weak linear increase observed for F/Fq > 0.30. We already 
mentioned that for such small F/Fq the positive signal near time zero is comparable in 
size to the demagnetization observed at later delays. It is therefore well possible that the 
obtained is influenced by the processes that are responsible for the positive A6/6q and 
we refrain from further conclusions based on Tm obtained for F/Fq < 0.30. 

We fitted the values for F/ Fq > 0.30 by a linear dependence. The result is shown in Fig. [3] 
by a line. We find that the demagnetization time increases by 67 ± 30 fs within F/Fq = 1 
with a nominal zero fluence limit of Tm = 690 ± 20 fs. 

For itinerant ferromagnets like Ni or Co an observed increase in Tm with fluence in com- 
bination with a finite value at zero fluence was explained by model descriptions based on 
Elliott- Yafet scatteringi^ and Stoner excitationsA^. The case of Gd has been discussed in 
Ref.— . It features two separate demagnetization times in agreement with the experimental 
observation—. These two time scales are attributed (i) to cooling of the hot electron sys- 
tem through e-ph coupling for the fast time scale, which is discussed in the present article, 
and (ii) to the demagnetization determined by the equilibrium spin-flip probability, which 
is weaker in Gd than for 3d transition metal ferromagnets and sets the slower timescale 
observed in Refs. 



This description predicts a fluence dependence of the ultrafast demagnetization of Gd 
that is determined by the increase of the transient lattice temperature Ti. We have cal- 
culated Te{t) and T\{t) for different fluence by the well known two-temperature model- 
in refined versions, which were published earlier— regarding its application to Gd(OOOl) 
see^. Considering a variation in fluence following the experimentally investigated range 
we find a pronounced variation of the time scale on which Te(t) and T\{t) equilibrate (not 
shown). This change in time scale is at a first glance more pronounced than the changes 
found for the demagnetization time in Fig. [31 However, a systematic analysis of the time 
scale at which Te{t) and T\{t) are changed is non-trivial since they vary in a non-exponential 
way. Therefore, we turn to a discussion of the excess energy e of the electronic system which 
is related to the electron temperature through e = 7Te(t)^, with 7 being the linear parameter 
in the temperature dependent specific heat of the electron system. Fig. H] shows the results 
obtained for e{t) with F/Fq varying between 0.2 and 1. The transient behavior can be de- 
scribed by a single exponential that represents energy transfer from the electron system to 
the lattice. The exponential time scales were determined with 100-800 fs and are plotted 
in the inset of Fig. |H The energy transfer time shifts to larger values for higher fluence and 
changes by three times in the investigated fluence range. Consequently the energy content 
of the lattice changes with these varying time constants. 

The comparison between time- dependent demagnetization curves and the time evolution 
of excess energy emphasizes two points which lead to the conclusion that the explanation of 
the ultrafast demagnetization in Gd solely through the time- dependent electron or lattice 
temperature is incomplete. At first we argue that the range of the energy transfer times r^-ph 
found in the investigated fluence range is with 0.2-0.6 ps considerably broader than the range 
observed for which is about 0.65-0.8 ps (see Figs. EE])- Second, and maybe more general, 
the demagnetization can be described by a simple exponential time dependence similar to 
the transient energy density. In contrast, the transient electron and lattice temperatures 
follow a more complicated evolution. Moreover, spin fluctuations which present a sizeable 
energy content of the system in particular in excited magnetic systems^^, might be essential 
to take into account. 
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FIG. 4: Main panel: Time-dependent excess energy density for different relative fluence, which 
was calculated by the two-temperature model. The inset shows the characteristic times Tq— ph of 
energy transfer from the electronic system to the lattice mediated by electron-phonon scattering 
which was determined by single exponential decay to the excess energy density in time interval 
from 100 to 800 fs. 

C. Discussion 

Our pump-probe analysis of fs laser-induced magnetization dynamics in Gd(OOOl) re- 
vealed two qualitatively different mechanisms, which both change transiently the magnetic 
moment per atom. The one mechanism is transfer of angular momentum from the spin 
system to another degree of freedom of the sample, which is finally the lattice. In this case 
the magnetic moment of the sample as a whole is reduced. Such mechanisms are discussed 
in the literaturei^ii^ to offer consistent descriptions of ultrafast demagnetization as observed 
in several time-resolved magneto-optical studies, see e.g. Ref. and x-ray magnetic-circular 
dichroism (XMCD) experiments^ for Ni. The observed demagnetization time of about 700 fs 
determined for Gd(OOOl) can be explained by relaxation of the optically excited electrons 
by interactions with phonons or spin waves which subsequently interact with the lattice. 
On a more microscopic level the demagnetization process in Ni and Gd must be different, 
because in a transition metal the magnetic moment is formed in the same band which is also 
optically excited. In a lanthanide like Gd the dominant part of the magnetic moment is local- 
ized near the ion core and resides in the 4/ level, which is not primarily optically excited^S. 
To achieve significant demagnetization the excess energy generated by the laser excitation 
must be transfered to the 4/ electrons. While details of this process will be discussed in a 
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forthcoming publication, we note here that optically excited 5d electrons are coupled to the 
localized 4/ electrons through intraatomic exchange interaction. Ultrafast changes of the 
magnetization, which are investigated by MOKE, can be expected to be dominated by the 
5d contribution to the magnetic moment. The femtosecond XMCD study performed at the 
Gd M5 edge probes the magnetic moment of the 4/ levels directly~. A comparison of the 
time scales of demagnetization obtained in both these experiments suggests that the 4/ and 
5d contribution to the magnetic moment are essentially strongly coupled on ultrafast time 
scales accessible here and demagnetize as a total magnetic moment. 

The second mechanism investigated here is spin transfer which we probe by means of 
the exchange-split spin polarized 5^22 surface state of Gd(OOOl). Here the magnetic mo- 
ment integrated over the sample remains constant, but is redistributed between different 
parts of the sample, in our case among surface and bulk electronic states. We explain this 
process to originate from the optical transitions that couple the exchange-split surface and 
bulk states. In the literature a similar scenario has been discussed recently as a potential 
explanation for bulk demagnetization in NiM. In this theoretical study the authors consider 
the optically excited transport of hot carriers in metallic samples which redistribute the 
magnetic moment from the ferromagnetic layer into a para- or diamagnetic substrate. Such 
transport effects of hot carriers are known since early investigations of femtosecond electron 
dynamics in metallic layera^^, however, to what fraction they are indeed responsible for the 
ultrafast demagnetization remains an interesting question which is to be clarified in future 
investigations. 



IV. CONCLUSIONS 



By combining ultrafast bulk and surface sensitive magneto-optical techniques we have 
analyzed the femtosecond laser-induced demagnetization of epitaxial Gd(OOOl) films. In the 
bulk the demagnetization occurs with a characteristic time of about 0.7 ps. Considering 
that the surface sensitive signal changes within the laser pulse duration of 35 fs and taking 
resonant optical transitions between valence electronic surface and bulk states into account, 
we attribute this ultrafast change to transfer of spin-polarized charge carriers between surface 
and bulk states. Variation of the pump fiuence up to 1 mJ/cm^ shows a weak increase in 
the bulk demagnetization time of about 70 fs/(mJ/cm^). A comparison with the fiuence 
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dependence expected from the energy transfer among electrons and the lattice revealed that 
albeit there is qualitative agreement in the increasing trend of demagnetization time such 
description remains incomplete due to pronounced quantitative variations. 
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